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ABSTRACT 

We report the discovery of a new Wolf-Rayet star in Aquila via detection of its cir- 
cumstellar nebula (reminiscent of ring nebulae associated with late WN stars) using 
the Spitzer Space Telescope archival data. Our spectroscopic follow-up of the central 
point source associated with the nebula showed that it is a WN7h star (we named it 
WR121b). We analyzed the spectrum of WR 121b by using the Potsdam Wolf-Rayet 
(PoWR) model atmospheres, obtaining a stellar temperature of ~ 50 kK. The stellar 
wind composition is dominated by helium with ~ 20 per cent of hydrogen. The stellar 
spectrum is highly reddened {Eb-v = 2.85 mag). Adopting an absolute magnitude 
of M v = —5.7, the star has a luminosity of \ogL/L Q — 5.75 and a mass-loss rate 



of 10 



-4.7 



M 



©yr" 



and resides at a distance of 6.3 kpc. We searched for a possible 



parent cluster of WR 121b and found that this star is located at ~ 1 degree from the 
young star cluster embedded in the giant Hil region W43 (containing a WN7+a/OB? 
star - WR121a). We also discovered a bow shock around the 09.5111 star ALS9956, 
located at ~ 0.5 degree from the cluster. We discuss the possibility that WR 121b and 
ALS 9956 arc runaway stars ejected from the cluster in W43. 

Key words: line: identification - circumstellar matter - stars: individual: ALS 9956 
- stars: Wolf-Rayet - open clusters and associations: individual: [BDC99] W43 cluster 



1 INTRODUCTION 

Wolf-Rayet (WR) stars are evolved massive stars possess- 
ing a strong stellar wind. Interaction of the WR wind with 
the material lost during the preceding evolutionary phase 
results in the origin of a compact circumstellar nebula (e.g. 
Chevalier & Imamura 1983; Robberto et al. 1993; Brighenti 
& D'Ercole 1995). The majority of WR circumstellar neb- 
ulae is associated with late WN-type (WNL) stars (Gvara- 
madze et al. 2009 and references therein), i.e. with stars that 
only recently entered the WR phase (e.g. Moffat, Drissen & 
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Robert 1989). These nebulae could be sources of mid- and 
far-infrared (IR) emission (e.g. Mathis et al. 1992; Hutse- 
mekers 1997; Marston et al. 1999) and can be detected with 
the modern IR surveys. Barniske et al. (2008) found circum- 
stellar emission from dust and gas around two extremely 
luminous WN stars near the Galactic center. 

The detection of compact nebulae reminiscent of nebu- 
lae associated with known WNL and Luminous Blue Vari- 
able (LBV) stars might suggest that their central stars are 
evolved massive stars as well and could be used for selection 
of candidate WNL, LBV or related stars for spectroscopic 
follow-ups. The importance of this channel for search for 
evolved massive stars was pointed out by Gvaramadze et 
al. (2009), who reported the discovery of a ring nebula in 
Cygnus (using the archival data from the Cygnus-X Spitzer 
Legacy Survey; Hora et al. 2009) and the results of spec- 
troscopic follow-up of its central star, showing that the star 
belongs to the WN8-9h subtype. Inspired by this discovery, 
we searched for more ring nebulae in the archival data of the 
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Figure 1. Spitzer MIPS 24 (im image of a new ring-like nebula in 
Aquila and its central star (WR 121b; located at the very centre of 
the ring). The signatures of a second incomplete concentric shell 
are indicated by three arrows. An arclike filament (surrounded by 
a rectangle) at ~ 3' to the north of the nebula is the south rim 
of the supernova remnant 3C 391. 



Spitzer Space Telescope and found dozens of objects (Gvara- 
madze, Kniazev & Fabrika 2009; cf. Carey et al. 2009) . 

In this paper, we present the results of study of one of 
the stars selected via detection of their circumstellar nebu- 
lae. Our spectroscopic follow-up of this star showed that it 
is a WR star of the WN7h subtype, which provides a fur- 
ther proof of the effectiveness of our method for revealing 
massive post-main-sequence stars. 



2 IR NEBULA AND ITS CENTRAL STAR 

The new nebula in Aquila was discovered in the archival data 
of the Spitzer Space Telescope, obtained with the Multiband 
Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004) 
within the framework of the 24 and 70 Micron Survey of 
the Inner Galactic Disk with MIPS (MIPSGAL; Carey et 
al. 2009). The MIPS 24 /xm image of the nebula (see Fig.[TJ» 
shows a perfect circular shell with a diameter of ~ 3.8 ar- 
cmin. One can see also signatures of a second incomplete 
concentric shell of diameter of ~ 4.6 arcmin and some struc- 
tures within the main shell; it is not clear, however, whether 
they are related to the nebula. 

The region containing the nebula was also covered by 
the Galactic Legacy Infrared Mid-Plane Survey Extraordi- 
naire (Benjamin et al. 2003) carried out with the Infrared 
Array Camera (IRAC; Fazio et al. 2004), the Multi- Array 
Galactic Plane Imaging Survey (Helfand et al. 2006), the 
SuperCOSMOS H-alpha Survey (SHS; Parker et al. 2005) 
and partially by Chandra observations of the (background) 
supernova remnant 3C 391 (located at ~ 3 arcmin to the 
north of the IR shell). None of these observations nor the 
MIPS 70 fj,m data show signatures of the nebula. 

Fig. [T] also shows a point source located exactly in the 
centre of the nebula, which suggests that this source might 
be associated with the nebula. The details of the source are 
summarized in Table[T] The coordinates and the J, H, K s 



Table 1. Details of the central star (WR121b) associated with 
the new nebula in Aquila 



RA(J2000) 


18 h 49 m 27?34 


Dec.(J2000) 


-01°04'20:'8 


l,b 


31.7515, -0.0480 


V (mag) 


17.15 ±0.03 


I (mag) 


13.60 ±0.03 


J (mag) 


10.94 ±0.03 


H (mag) 


10.07 ±0.03 


K s (mag) 


9.47 ±0.02 


[3.6] (mjy) 


77.8 ± 2.3 


[4.5] (mjy) 


70.6 ±2.1 


[5.8] (mjy) 


55.5 ± 1.7 


[8.0] (mjy) 


37.6 ± 1.1 


[24] (mjy) 


9.7 ±0.3 


[70] (mjy) 


< 3.1 



magnitudes are taken from the 2MASS All-Sky Catalog of 
Point Sources (Skrutskie et al. 2006). 

The photometric V and I magnitudes of the optical 
counterpart to the IR central source were determined on 
a CCD frame obtained with the SBIG CCD ST-10XME at- 
tached to the 40-cm Meade telescope of the Cerro Arma- 
zones Astronomical Observatory of the Northern Catholic 
University (Antofagasta, Chile) during our observations in 
2009 April. The accuracy of the measurements is ~ 0.03 mag 
in both filters. 

To measure the Spizer's IRAC and MIPS fluxes from 
the star, we performed its aperture photometry using the 
M0PEX/APEX source extraction package and applied an aper- 
ture correction obtained from nearby bright point sources. 
The estimated errors are ~ 2 — 3 per cent. For the flux at 
70 /im, we give a 3a upper limit since the star was not de- 
tected at this wavelength. 



3 SPECTROSCOPIC OBSERVATIONS AND 
DATA REDUCTION 

To determine the nature of the central star associated with 
the new nebula in Aquila, we obtained its spectrum with the 
TWIN spectrograph attached to the Cassegrain focus of the 
3.5 m telescope in the Observatory of Calar Alto (Spain) dur- 
ing director's discretionary time on 2009 May 5. The setup 
used for TWIN was the grating T08 in the first order for the 
blue arm (spectral range 3500-5600 A) and T04 in the first 
order for the red arm (spectral range 5300-7600 A) which 
provides with an inverse dispersion of 72 A/mm for both 
arms. We have used the CCD detectors SITe#22b_14 and 
SITe#20b_12 for the blue and red arms, respectively, with 
the 5500 A beam splitter. A slit of 240" x 2.1" was used for 
these spectral observations. The resulting FWHM spectral 
resolution measured on strong lines of the night sky and 
reference spectra was 3-3.5 A. The seeing during the ob- 
servations was stable, ~ 1.5 — 1.6 arcsec. The observations 
were carried out at relatively low airmasses (1.25—1.40) and 
the slit was aligned exactly to the parallactic angle to ex- 
clude any atmospheric dispersion effects. The total exposure 
of 3600 s was broken up into three 1200 s subexposures to 



New Wolf-Rayet star in Aquila 3 



■13.5 



■14 



14.5 



-15 - 



GO 

° -15.5 



■16 




4500 



5000 



5500 



6000 



6500 



7000 



Wavelength (A) 



Figure 2. The observed spectrum of the central star (WR 121b) associated with the new nebula in Aquila with the most prominent 
lines indicated. 



allow for removal of cosmic rays. Spectra of He-Ar compar- 
ison arcs were obtained to calibrate the wavelength scale. 
The spectrophotometric standard star BD+33° 2642 (Oke 
1990; Bohlin 1996) was observed at the beginning for flux 
calibration. 

The primary data reduction was done using the IRAF0 
package ccdred. The data for each CCD detector were 
trimmed, bias subtracted and flat corrected. We used the 
IRAF software tasks in the twodspec package to perform the 
wavelength calibration and to correct each frame for distor- 
tion and tilt. To derive the sensitivity curve, we fitted the 
observed spectral energy distribution of the standard star 
by a low-order polynomial. The final sensitivity curves have 
an internal precision better than 1 per cent over the whole 
optical range. One-dimensional (ID) spectra were then ex- 
tracted using the IRAF APALL task. All one-dimensional 
spectra obtained with the same setup were then averaged. 
Finally, the blue and red parts of the total spectrum were 
combined without additional factors. The resulting reduced 
spectrum is shown in Fig. [2] 

1 IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatory, which is oper- 
ated by the Association of Universities for Research in Astron- 
omy, In. (AURA) under cooperative agreement with the National 
Science Foundation (NSF). 



The main emission lines were measured applying the 
MIDAS programs described in detail in Kniazev et al. (2004, 
2008a). Their equivalent widths (EW), FWHM and helio- 
centric radial velocities are summarized in Table[2] 



4 A NEW WOLF-RAYET STAR - WR 121B 
4.1 Spectral type 

Fig. [2] shows that the spectrum of the central star of the 
new nebula in Aquila is dominated by strong emission lines 
of hydrogen, He I, Hell, Nni and Niv. A less pronounced 
emission line of C IV is also present in the spectrum. Many 
of the weaker lines show P Cygni profiles, while the strongest 
lines are entirely in emission. Numerous absorptions visible 
in the spectrum are either interstellar [e.g. diffuse interstellar 
bands (DIBs) at AA4428, 4726, 4762, 5780, 5797 and 6280] or 
(for A > 7200 A) telluric in origin. No forbidden lines can be 
seen in the spectrum. 

The prominent helium and nitrogen emission lines in- 
dicate that the star belongs to the nitrogen sequence of 
Wolf-Rayet stars (WN). To determine its WN subtype we 
apply the three-dimensional classification scheme by Smith 
et al. (1996). The measured line ratios Hell A5412/Hel 
A5876, Civ A5808/Hen A5412 and Civ A5808/Hei A5876 
(~ 0.9, 0.2 and 0.2, respectively) suggest that the star be- 
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Table 2. Equivalent widths (EW), FWHM and heliocentric radial 
velocities of main emission lines. 



A (A) Ion 



EW(A) 
[A] 



FWHM(A) 

[A] 



[km s" 1 ] 



4541 


He II 


5.3±2.3 


12.4±0.8 


77±11 


4640 


N m 


27.1±6.6 


22.7±0.9 


131±11 


4686 


He II 


33.8±4.0 


13.9±0.3 


124±6 


4861 


He ll+H/3 


10.2±2.8 


12.3±0.6 


143±8 


5412 


He II 


8.9±1.0 


10.7±0.2 


142±5 


5808 


C iv 


2.5±0.6 


20.0±2.2 




5876 


He I 


11.2±1.0 


14.0±1.0 


164±10 


6234 


He II 


0.9±0.2 


7.9±0.4 


148±6 


6311 


He II 


0.4±0.1 


4.2±1.3 


132±11 


6406 


He II 


1.6±0.3 


10.0±0.6 


130±8 


6563 


He ii+Hq 


63.2±0.7 


21.3±0.1 


133±4 


6678 


He i+Hc II 


14.9±0.7 


27.3±0.5 




6891 


He II 


2.9±0.5 


10.6±1.8 


186±12 


7065 


He I 


11.5±0.8 


29.8±2.0 




7111 


N iv 


14.2±0.5 


14.4±0.9 


140±9 


7123 


N iv 


14.3±0.5 


12.8±0.7 


120±8 


7178 


He II 


3.8±0.3 


9.6±0.6 


142±10 


7281 


He I 


0.8±0.1 


6.3±1.0 


129±12 



longs to the WN7 subtype. The same classification follows 
from the position of WR 121b in the plot of He II A5412 EW 
versus He n A4686 FWHM [see Fig. 13 and Fig. 15 of Smith 
et al. (1996) and Table[2]. 

The hydrogen Balmer lines coincide with the lines of 
the He II Pickering series 4 — n for even principle quantum 
numbers n. The quantitative analysis (see below) will reveal 
that the Hell blends with Ha and H/3 are stronger than 
predicted by a hydrogen-free model that fits the unblended 
Hell 7-4 line at 5412 A. Hence there is a detectable contri- 
bution from hydrogen to these line blends. This conclusion 
can be justified by using the empirical 'hydrogen criteria' 
for WN stars introduced by Smith et al. (1996). For EWs of 
lines at 4861, 4541 and 5412 A given in Table[2] one finds a 
ratio EW(4861)/[^EW(4541)EW(5412) - 1] = 1.7, which 
is much larger than 0.5, the figure defining the presence of 
hydrogen. Following the three-dimensional classification for 
WN stars by Smith et al. (1996), we therefore add a suffix 
'h' to the spectral subtype, which then becomes WN7h. 

We name our program star WR 121b, in accordance 
with the numbering system of the Vllth Catalogue of Galac- 
tic Wolf-Rayet Stars by van der Hucht (2001). 



4.2 Nebular emission 

The red part of 2D spectrum shows a complex of emis- 
sion lines originating in the region covered by the slit. 
The intense emission lines include, in particular, Ha, 
[Nil] A6584, and [Sn] AA6716, 6731. The ratios of the 
line fluxes (log([N n]/Ha)=-0.73, log([Sn]/Ha)=-0.58, 
[S ll](6717/6731)=1.3) clearly show that the lines originate 
in a hot tenuous medium, which is most likely an extended 
Hll region [see, e.g., the diagnostics of Hn regions by Kni- 
azev et al. (2008b)]. The intensities of all these lines are very 
uniform along the slit. The Ha radial velocity [measured in 
the way described by Zasov et. al. (2000)] is also constant 
(36 ± 4 kms -1 ) along the slit. We conclude therefore that 



the nebular line emission originates from an Hll region un- 
related to WR 121b and its ring nebula. 

This conclusion is supported by an estimate of the in- 
terstellar extinction from the Balmer decrement. We found 
Ay = 4.6 mag (constant along the slit), which is much less 
than Ay = 8.8 mag derived from the stellar spectrum (see 
Sectionl4~4l) . 



4.3 Spectral analysis and stellar parameters 

To analyze the stellar spectrum and to derive the funda- 
mental parameters of WR 121b, we use the Potsdam Wolf- 
Rayet (PoWR) models for expanding stellar atmospheres. 
These models acount for complex model atoms including 
iron-line blanketing in non-LTE (for a detailed description 
see Hamann & Grafener 2004) . For abundances of trace ele- 
ments we adopt mass fractions which are typical for Galactic 
WN stars - N: 0.015, C: 0.0001, Fe: 0.0014 per cent (Hamann 
& Grafener 2004). 

The main parameters of a WR-type atmospheres are 
the stellar temperature, T«, and the so-called transformed 
radius, R t . The stellar temperature T* denotes the effective 
temperature related to the radius R t , i.e. L = aTfATvRi, 
where a is the Stefan-Boltzmann constant and i?* is by def- 
inition at a Rosseland optical depth of 20. Rt is related to 
the mass-loss rate M and defined by 



Rt — R* 



2500 kms- 



10 -4 M yr- 



2/3 



(1) 



where D is the clumping contrast, and v x is the terminal 
velocity of the wind. 

These basic stellar parameters T* and Rt are deter- 
mined from fitting the lines in the normalized spectrum (see 
Fig-El • The normalization is in fact achieved by dividing the 
absolutely calibrated observed spectrum by the theoretical 
continuum, which makes the total procedure described in 
the following an iterative process. 

A first orientation about the proper choice of these pa- 
rameters can be obtained by comparing the observed spec- 
trum to the published grids of WN models (Hamann & 
Grafener 2004) . With respect to the stellar temperature, we 
select the model which reproduces best the balance between 
the lines from He I versus He II. At the same time, Rt is prop- 
erly chosen which influences the strength of the emission 
lines in general. Conveniently, the grid model WNL07-11 [T* 
= 50 kK, log(7? t /^©) = 1-0] provides a satisfactory fit to the 
observed line spectrum (see Fig.|3}. Note that according to 
its location in the transformed radius-temperature diagram 
(Hamann & Grafener 2004), WR121b belongs to the WN7 
subtype, consistent with our classification in Section l4~T1 

The terminal wind velocity of = 1000 kms -1 , for 
which the grid had been calculated, nicely reproduces the 
widths of the line profiles and thus obviously adequate for 
our star. 

The clumping contrast D [cf. equation (1)] has influence 
on the electron-scattering wings of strong emission lines. 
Hamann & Koesterke (1998) established D — 4 as a typ- 
ical value for WN stars. The model grid from Hamann & 
Grafener (2004) which we use here is calculated with D — 4, 
and therefore we keep this choice for practical reasons. Re- 
cent discussions (see Hamann et al. 2008) favour stronger 
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Figure 3. Observed optical spectrum (blue/solid line) of WR121b, compared with the best-fitting model (red/dotted line) with the 
parameters as given in Table|3] For normalization the absolutely calibrated observation was divided by the reddened model continuum, 
and after slightly aligned further "by eye" . 



clumping up to D — 10. Data and fit quality do not allow 
an individual estimate of D for WR 121b to this precision, 
but values smaller than 4 can be ruled out e.g. from the red 
wing of Ha/He II 6-4. Note that the effect of D is simply a 
scaling of the empirically derived mass-loss rate proportional 
to D~ 1/2 . 

The WNL model grid is calculated with a hydrogen 
abundance (Xh) of 20 per cent by mass. For a closer check 
we compare the fit with models of zero hydrogen and with 
40 per cent hydrogen. Fig.[4]shows the lines of the He II Pick- 
ering series (8-4: 4862 A; 7-4: 5412 A; 6-4: 6562 A). Those 
profiles that are blended with H/3 (left frame) or Ha (right 
frame) become much stronger with growing hydrogen abun- 
dance, while the unblended Hell line at 5412 A (middle 
frame) becomes slightly weaker because of the compensat- 
ing decrease of helium. The comparison with the observed 
profiles (blue/thick solid lines in Fig. [4} reveals that the 
hydrogen-free model predicts a too small line at the Ha 
wavelength, while with 40 per cent hydrogen Ha and H/3 
are both stronger in the model than observed. Hence the 



hydrogen abundance of 20 per cent is a good estimate, and 
we can consider the grid model as our best fit. 



4.4 Luminosity and distance of WR 121b 

The spectral analysis alone cannot tell the absolute dimen- 
sions of the star. These are related to the distance of the 
object, which is a priori unknown for WR121b. 

However, a couple of Galactic WN stars can be assigned 
to open clusters or associations for which the distance can 
be independently determined. On this basis, correlations be- 
tween the absolute magnitude and the spectral subtype have 
been established. 

Hamann et al. (2006) have shown that stars with WN7 
subtype come in two flavours. Some WN7 stars belong to 
the group of extremely bright and luminous WNL stars that 
come from very massive progenitors. Their absolute magni- 
tude is about M v — —7.2 mag (v magnitudes refer to the 
monochromatic flux at 5160 A), and their winds typically 
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Figure 4. Test for the hydrogen abundance in WR 121b. The observed profiles (blue/thick solid line) are compared with models (red/thin 
lines) for different hydrogen abundance, but otherwise the same parameters as of our best-fit model (see Tablel3l . The tested hydrogen 
mass fractions are zero (dashed lines), 20 per cent (thin solid lines) and 40 per cent (dotted lines), respectively. 



show a high hydrogen abundance of about 50 per cent by 
mass. 

However, stars of the spectral subtype WN7 can also be- 
long to an extension of the "early" subtypes (WNE), which 
are less luminous. According to the subtype calibration by 
Hamann et al. (2006), a WN7 star of this kind would have an 
absolute magnitude of M v — — 5.7 mag (with an uncertainty 
of about 0.5 mag). Most of the WNE stars are hydrogen-free, 
but some of them show a small hydrogen mass fraction be- 
low 25 per cent. According to the hydrogen content of about 
20 per cent determined above, WR 121b rather belongs to 
this group of less bright stars. 

Therefore we now scale our best-fit grid model to an ab- 
solute brightness of M v = —5.7 mag. Remember that, in fair 
approximation, WR models of the same "transformed ra- 
dius" Rt [see equation (JTJ] can be scaled to different absolute 
sizes. The scaled model has a luminosity of log L/Lq = 5.75 
(cf. Table|3J). The corresponding stellar radius, mass-loss rate 
and the number of hydrogen-ionizing photons ($i) are in- 
cluded in the same table. Note that the error estimates given 
in Table[3] are rough and based on the experience that the 
line fit is accurate to one cell of the model grid, i.e. to 0.05 
in logT and to 0.1 in logi? t . The main uncertainty in the 
luminosity and the distance comes from the subtype calibra- 
tion of the absolute visual magnitude. The given mass-loss 
rate refers to the adopted clumping parameter of D = 4. 
The uncertainty in the mass-loss rate is mainly due to the 
error margin of the luminosity, since the empirical M scales 
proportional to L 2//3 . 

The spectral energy distribution (SED) of the scaled 
model is now fitted to the photometric observations (Fig. [5}. 
Two parameters can be adjusted for the fit, the distance d 
and the reddening Eb-v- The reddening law we adopt from 
Seaton (1979) in the optical and from Moneti et al. (2001) 
in the IR. A perfect fit to the whole SED is achieved with 
Eb-v = 2.85 mag and d — 6.3 kpc. The strong interstellar 
absorption also shows up in the observed spectrum by pro- 



Table 3. Stellar parameters for WR 121b 



Spectral type 


WN7h 


T» [kK] 


50±5 


\ gRt[R & \ 
Voo [kms -1 ] 


10±4 


1000±200 


X H 


0.2±0.1 


log L [Lq] 


5.75±0.30 


R* [Rq] 


10.0±2.0 


logMIMeyr- 1 ] 


-4.7±0.3 


log*; [s- 1 ] 


49.6±0.3 


E B -v [mag] 


2.85±0.05 


d [kpc] 


e Q + 2.2 

D -°-i.e 



nounced DIBs (see Fig. [2}. The high reddening implies an 
extinction in the visual of Ay = 8.8 mag and a severe de- 
crease of the flux at shorter wavelengths; for A < 4400 A our 
observed optical spectrum drowns in the noise. 



5 POSSIBLE BIRTH CLUSTER OF WR121B 

WR 121b, like the majority of WR stars, is located outside of 
known star clusters. Since most of massive stars originate in 
a clustered mode of star formation (e.g. Zinnecker & Yorke 
2007), it is natural to assume that WR121b (or its progeni- 
tor star) was ejected from the parent cluster either through 
a dynamical process in the core of the cluster (Poveda et al. 
1967; Gies & Bolton 1986) or due to a supernova explosion 
in a close massive binary system (Blaauw 1961; Stone 1991). 
The first process could start to operate at the very beginning 
of cluster evolution, while the second one only several Myr 
after the cluster formation, when the most massive stars in 
binaries (and respectively in the cluster) end their lives as 
supernovae. In both cases, the ejected WR star should not 
be far from its parent cluster since the progenitors of WR 
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Figure 5. Observed flux distribution of WR 121b (blue/noisy) in 
absolute units, including the calibrated spectrum and the photo- 
metric measurements compiled in Tablc[T] compared to the emer- 
gent flux of the model continuum (red/smooth line), in the op- 
tical also shown with lines. The model flux has been reddened 
and scaled to the distance according to the parameters given in 
Tablc[3] 

stars are very massive (> 25 Mq) and therefore short-lived 
(< 5 Myr) stars. 

Proceeding from the aforementioned, we searched for 
known young star clusters around WR 121b using the SIM- 
BAE0 and WEBDA0 databases and found that this star 
is located at ~ 1 degree from a cluster deeply embedded 
(Av ^ 30 mag) in the giant Hn region W43 (see Fig.|6j). 
This cluster (named in the SIMBAD database as ' [BDC99] 
W43 cluster') contains at least three evolved massive stars 
(Blum, Damineli & Conti 1999), two of which are O-type 
giants or supergiants and the third one is a WN7+a/OB? 
star (WR121a). 

The large ionizing flux of the W43 star forming region 
(~ 10 51 Lyman continuum photons s , Smith, Biermann 
& Mezger 1978) is comparable with that of the very mas- 
sive star cluster NGC 3603, which suggests that the central 
cluster in W43 contains a large number of yet undetected 
massive stars and therefore should be effective in producing 
runaway stars (e.g. Gvaramadze 2009; Gvaramadze, Gua- 
landris & Portegies Zwart 2008, 2009). On the other hand, 
the presence of the WNL star in the cluster implies that the 
age of the cluster is comparable to the evolutionary age of 
WR 121b. It is plausible, therefore, to identify the origin of 
WR 121b (or its progenitor star) with the central cluster of 
W43. 

This conjecture is also supported from the distance to 
WR 121b of 6.3± \% kpc which we have obtained in Section 
14.41 For W43 a distance of about 7 kpc was derived from H I 
radio measurements and a Galactic rotation model (Wilson 
et al. 1970), which for the modern value of the distance to 
the Galactic Centre reduces to about 5.6 kpc (Blum et al. 
1999). Blum et al. (1999) obtained 5.7 kpc from two O stars 
observed in W43 when assuming that they are giants. Thus, 
within the error margins the distances are consistent. 

With a luminosity of log L/Lq ~ 5.75, the current and 
the initial masses of WR121b are ~ 20 and 40 Mq, respec- 
tively, and its evolutionary age is ~ 4 — 5 Myr. The latter 

2 http:/ /simbad. u-strasbg.fr/simbad/ 

3 http://www.univie.ac.at/webda/ 
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Figure 6. A 2° X 2° 8.28 /jm image of the Hn region W43 and 
its environments obtained in the Midcourse Space Experiment 
Galactic plane survey (Price et al. 2001), with the position of the 
star cluster embedded in W43 marked by a white circle and the 
positions of WR 121b and ALS 9956 indicated by diamond points. 



See text for details. 
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Figure 7. Left: Spitzer MIPS 24 /an image of ALS 9956 (marked 
by a circle) and its bow shock. Right: SHS image of the same field. 
Orientation of these images is the same as for Fig. [6] 



estimate and the angular separation between WR 121b and 
W43 (~ 1 degree, which for d — 6.3 kpc corresponds to 
the linear separation of ~ 110 pc) implies that the mini- 
mum allowed peculiar (transverse) velocity of WR 121b is 
~ 25 kms -1 . We predict, therefore, the 'observed' (i.e. not 
corrected for the Galactic rotation and the solar motion) 
proper motion of WR 121b of [i a ~ — 1.6 mas yr -1 and 
Us — — 3.7masyr _1 (if the star was ejected ~ 4 — 5 Myr 
ago) or fi a ~ — 0.3 mas yr -1 and {is ~ — 1.2 mas yr -1 (if the 
the ejection occurred ~ 1 Myr ago). [These estimates were 
derived using the Galactic constants Rq — 8 kpc and G = 
200 kms -1 (e.g. Reid 1993; Avedisova 2005) and the solar 
peculiar motion (Uq,Vq,W q ) = (10.00,5.25,7.17) kms" 1 
(Dehnen & Binney 1998).] Proper motion measurements for 
WR 121b are necessary to prove the association between this 
star and the star forming region W43 and thereby to con- 
strain some of its inferred fundamental parameters. 
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To find the evidence that the [BDC99] W43 cluster loses 
its massive stars, we searched for bow shock-like structures 
in its environments using the MIPSGAL survey. The visual 
inspection of MIPS 24 fim images results in the discovery 
of a bow shock produced by the 09.5111 star ALS 9956 (see 
Fig.0, located at ~ 0.5 degree from the cluster (see Fig.|SJ). 
The bow shock has an obvious optical counterpart in the 
SHS. The orientation of the symmetry axis of the bow shock 
and the evolutionary age of the star of ~ 4 — 5 Myr are 
consistent with a possibility that ALS 9956 was ejected from 
the cluster in W43 (cf. Gvaramadze & Bomans 2008). 

Further support to this possibility comes from compar- 
ison of the (photometric) distance to ALS 9956 with the dis- 
tance to the cluster. ALS 9956 has the B and V magnitudes 
of, respectively, 11.83 and 11.25 (Kharchenko 2001). Using 
the extinction law from Rieke & Lebofsky (1985) and the 
synthetic photometry of Galactic O stars by Martins & Plez 
(2006), we derived extinction towards ALS 9956 of Ay ~ 2.6 
mag and estimated the distance to this star to be 6.2 kpc, 
i.e. similar to the distance to the cluster. 

Using the VizieR catalogue access tooQ, we searched 
for proper motion measurements for ALS 9956 and found 
that all existing measurements are insignificant (i.e. the 
measurement uncertainties are comparable to the measure- 
ments themselves) and, therefore, cannot be used to prove 
unambiguously our hypothesis that the star was ejected from 
W43. Still, taken at face value, they suggest that ALS 9956 
is a runaway star moving away from the Galactic plane to- 
wards growing Galactic latitude (which is consistent with 
the position of the star above the Galactic plane), while the 
direction of stellar motion inferred from the symmetry of the 
bow shock (see, e.g., Gvaramadze & Bomans 2008) is well 
within the stellar proper motion error cone. 

We note that the visual extinction towards ALS 9956 
and WR 121b is much smaller than that towards the 
[BDC99] W43 cluster. If both stars were indeed ejected from 
the cluster, then we conclude that some OB stars can be 
detected in optical wavelengths only because they are run- 
away stars, while their cousins residing in the parent clusters 
might still remain totally obscured. This may explain why 
some field O-type stars show no apparent association with 
star clusters [cf. de Wit et al. (2005); Schilbach & Roser 
(2008); Gvaramadze & Bomans (2008)]. 



6 CONCLUSION 

We have discovered a new WNL (WN7h) star via detec- 
tion of a circular IR nebula (typical of evolved massive 
stars) and spectroscopic follow-up of its central star. This 
discovery provides further evidence that the circumstellar 
nebulae associated with WR stars are inherent almost ex- 
clusively to WR stars of the WNL subclass. It also shows 
that the IR imaging could be an effective tool for reveal- 
ing WNL and related evolved massive stars. The new WNL 
star is one of many dozens of (candidate) evolved massive 
stars revealed via their association with compact IR nebu- 
lae, detected in the MIPS 24 /im data from the Spitzer Space 
Telescope archive. The spectroscopic study of these stars is 

4 http:/ /webviz.u-strasbg.fr/viz-bin/VizicR 



currently underway and its results will be presented in the 
forthcoming papers. 
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